PRESSURE REGULARITY CRITERION FOR THE THREE DIMENSIONAL 
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Abstract. In this paper we consider the three— dimensional Navier-Stokes equations in an infinite channel. We 
provide a suHicient condition, in terms of 9zp, where p is the pressure, for the global existence of the strong 
solutions to the three— dimensional Navier-Stokes equations. 
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1. Introduction 

The question of global regularity for the 31? Navier-Stokes equations is a major open problem in applied 
analysis. Over the years there has been an intensive work by many authors attacking this problem (see, e.g., 
[11], [H], [14], [22], [25], [27], [28], [29], [37], [40], [41], [42] and references therein). It is well known that the 2D 
Navier-Stokes equations have a unique weak and strong solutions which exist globally in time (cf., for example, 
[12j . [30], [H]). In the 3D case, the weak solutions are known to exist globally in time. But, the uniqueness, 
regularity, and continuous dependence on initial data for weak solutions are still open problems. Furthermore, 
strong solutions in the 31? case are known to exist for a short interval of time whose length depends on the 
physical data of the initial-boundary value problem. Moreover this strong solution is known to be unique (cf. , for 
example, [12], [37], [40]). The subtle difference between the 2D and 3D incompressible Navier-Stokes equations 
manifests itself in a clear way in the vorticity formulation of these equations. Let u denote the velocity field and 
a; = V X u the vorticity. The equations that govern the evolution of the vorticity are given by 

duj 

— vl\uj + (m • V)uj — uj ■ Vu = 0, 
V-u = 0. 

In this formulation the main obstacle for proving the global regularity is the vorticity stretching term u ■ Vu. 
This term is identically equal to zero in the 2D case. Nonetheless, there are some results regarding the global 
regularity for the 3-D Navier-Stokes equations under special symmetries and for which the vorticity stretching 
term is non-trivial. For example, the case of 3D axi-symmetric flows in domains of revolution (with a positive 
distance from the z— axis) [22] and [23] and helical flows [26]. It is worth mentioning, however, that the questions 
of the global regularity for 3-D Euler equations (i.e., when the viscosity v = Q) for axi-symmetric and helical flows 
with nontrivial vorticity stretching term are still open. 

In addition, there are few known results regarding the global regularity for special types of initial data. For 
instance, for small initial data (cf. e.g., [T^], [30] )• Fujita and Kato IE proved the global well-posedness 
for small H'^/'^ initial data, and later Kato [20] proved the same result for small initial data (see also |17j). 
An interesting global existence result for large, but very "oscillatory", initial data was also proved in [3], [S], 
[5], [30] (see also references therein). The latter case can be regarded, roughly speaking, as small initial data 
in H^/^ and hence falls as a special case of results in [TSl. The proofs of the above result relies heavily on the 
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viscosity mechanism. Roughly speaking, the viscosity acts very strongly in dissipating the solution starting from 
such oscillatory initial data so that in a very short time the solution becomes very small. As a result from that 
moment on a small initial data argument applies to establish the global existence. Here again, it is worth stressing 
that nothing is known about global regularity for the 31? Euler equations for this kind of oscillatory initial data. It 
is also worth mentioning the result of ^34. about global existence of the strong solutions to the three-dimensional 
Navier-Stokes equations for large initial data in thin domains, where the upper bound on the size of the initial 
data depends inversely on the thickness of the domain. 

Taking a physical point view, we define the dimensionless Reynolds number Re — where Uq represents a 
typical magnitude of velocity (e.g. the size of the initial data) and £ a typical length scale of the domain (e.g. 
its thickness). As a result one can view, roughly speaking, the above mentioned global regularity results to hold 
under the assumption of small enough Reynolds number. 

The other direction concerning the question of global regularity for the 3D Navier-Stokes equations is to 
provide sufficient conditions for the global regularity. For example we refer the reader to the pioneer work of 
Prodi 33J and of Serrin (see also the survey paper of Ladyzhenskaya [53] and references therein). Most 
recently, there has been some progress along these lines which states that a strong solution u exists on the time 
interval [0, T] as long as 

u e i"([0,T],L^(17)), with - + 4 = 1, for /3>3, (1) 

a p 

(see, for example, [2], [3], [T7], [H], EH], [Mj, [SH], and references therein). Moreover, there are some sufficient 
regularity conditions only on one component of the velocity field of the 3-D NSE on the whole space or under 
periodic boundary conditions (cf. e.g., [19j, [21], [32], [43]). In [7] we introduced a sufficient regularity condition 
on one component of the velocity field of the 3Z) Navier-Stokes equations under Dirichlet boundary conditions 
(see also, [?T]). 

The Navier-Stokes equations of viscous incompressible fluid in an infinite channel x (0, 1) C read: 
dv 

— - I^AhV - VV^:, + (v ■ Vh)v + WVy^ + VhP = /, (2) 

— - J^A^w - vw^^ + V ■ VhW + wwz ^ g, (3) 

Vh-v + w, = 0, (4) 

where v = {vi,V2), the horizontal velocity field components, w, the vertical velocity component, and p, the 
pressure are unknowns. (,f,g), the forcing term, and v > 0, the viscosity are given. We set \7h = {dx,dy) to be 
the horizontal gradient operator and Ah — + dy is the horizontal Laplacian. We denote by: 

r„-{(a:,y,l)eM^}, (5) 
r6 = {(a;,y,0)eM3}, (6) 

the physical solid boundaries of the channel M."^ x (0, 1). We equip the system ([2])-((4]) with the following no-normal 
flow and stress-free boundary conditions on the physical boundaries of the channel, namely, 

= 0. ui = 0; on r„ and Tb- (7) 

oz 

Furthermore, we assume that (u, w) is periodic with period 1 in both horizontal directions. That is, 

v{x + 1, y, z) = u(x, y + l,z)^ v{x, y, z), w{x + 1, y, z) = w{x, ?; + 1, z) = w{x, y, z). (8) 
Because of the horizontal periodic boundary conditions we consider here 

fi = (0,1)3 
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as our basic domain of the flow, and we denote by 

M = (0,1)2. 

In addition, we supply the system with the initial condition: 

v{x,y,z,Q) ^ vo{x,y,z), (9) 
w{x,y,z,Q) wo{x,y,z). (10) 

Let us denote by L'^iyt) and H™{il) the usual L'^— Lebesgue and Sobolev spaces, respectively (cf. [T). We 
denote by 

\\(t)\\^ = (^J dxdydzy , for every </) G ^^(f}). (11) 

Let 

Vi = |0e(C-(R2x[O,l]))^:|^[_^ = O;|^ 
V2 = {ibeC^iR^ X [0,1]) : ^\^^a^O;^\^^,=0;i^ix + l,y,z) = i,ix,y + l,z) = ^ix,y,z)}, 

V = {(0,^) Vi, V2, V,,-0 + ^, = 0}. 

We denote by H and V be the closure spaces of V in L'^{il.) under the L^— topology, and in iJ^(ri) under the 
i/^— topology, respectively. We say {v,w) is a Leray-Hopf weak solution to the system (P1)- PU)) if {v,w) satisfies 

(1) iv,w) G C{[0,T], H) n L'^{[0,T],V), and {dtv,dtw) G L^{[0,T],V'), where V is the dual space ofV; 

(2) the weak formulation: 



= 0; 0(x + 1, y, z) = 4i{x, y + l,z) = 4>{x, y,z)} , 

z = l 



{v ■ (/) + wip) dxdydz — / (f (io) • (f'ito) + 'w(io) '/'(^o)) dxdydz 
Jn 

iv,w) ■ {{cj)t,ipt) + i'{^h'f> + 4>zz, ^hi^ + i'zz)) dxdydz ds 



[{v ■Vh){4>,'^) ■ iv,w) +w{4>^,i;^) ■ {v,w)] dxdydz + / / [{f , g) ■ {v,w)] dxdydz, 
to J^i Jto Jn 

for every (0, ^p) G V, and almost every t, to G [0, T]; 
(3) the energy inequality: 

lci(||H|i+||z«|||) ^^^||^^^||2^||^^^||2^|j^^||2^||^^|j2] < r (^f^g).^y^^) dxdydz. (12) 



2 dt 

Moreover, a weak solution is called strong solution of (Pl)- pT7)) on [0, T] if, in addition, it satisfies 

(v, w) G C{[0, T], V) n L\[0, T],H\n)). 

Remark 1. Notice that one can extend v and w to be periodic, vertically, by, first, setting 

v{x,y,z,t) ^ v{x,y,-z,t), zg(-1,0) 
w{x,y,z,t) = -w{x,y,-z,t), z G (-1,0), 

and then, setting 

v{x,y, z + 2,t) = v{x,y, z,t) w{x,y,z + 2,t) — w(x,y,z,t) Vz G M. 
Similarly, p, /, and vq can be extended to be even periodic functions in the z— vertical with period 2, and g and 
wq are odd periodic functions in the z— vertical with period 2. In this way, the system (|2|l- (fT0|) can be treated 
as the 3D Navier-Stokes equations with periodic boundary condition. In other words the system (P|)- pil)) is 
a special case of the ZD Navier-Stokes equations on M'^ with periodic boundary condition. Then, by standard 
procedure for the 2>D Navier-Stokes equations with periodic boundary condition (see, e.g., [12], [H], [IS], [M] . 
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|37j . [40] . [41] ) one can show that there exists, global in time, a Leray-Hopf weak solution to the system ([2|)- (fT0|) 
if {vo,wo) £ H. Furthermore, one can show the short time existence of the strong solution if (vq,wo) G V. 



In this paper, we provide sufficient conditions on the pressure which guarantee the global existence of the 
strong solution to the 3D Navier-Stokes equations in infinite channel subject to the boundary conditions IT])-®. 
Several authors (see, for example, [3], [10], [13], [35], [44], [J^) have studied the question of the global regularity 
of the 3D Navier-Stokes equation by providing sufficient conditions on the pressure. Specifically, the authors of 
[35j have shown that if |up + 2p is bounded from below or from above then the weak solution is strong solution. 
The authors of [10 have shown that if p satisfies 

L"(0,r;L'='(M3)) with - + 4<2, for /? > ^, (13) 

a p 2 

then the weak solution is strong solution. This result has been improved in [5] by assuming 

p e L"(0, T; L^(R3)) with - + I = 2, for /3 > % (14) 

a p 2 

In [IS] the author has established the global regularity by assuming that 

V,,p, eL"(0,r;L'3(M3)) ^iti^ - + |<3, for a > ^, /3 > 1. (15) 

a p 3 

Here, we will show the existence of the strong solutions of the system (j^- pHj) on interval [0, T] provided the 
vertical derivative of the pressure satisfies 

p, £ L°'{[0,T];L^{Q)) with a > 3, and /3 > 2. (16) 

Let us observe that the quantity that appears in the Navier-Stokes equations is {VhP,Pz) and not p itself. 
Therefore, the conditions from [45j and our conditions seem to be natural. Furthermore, it is worth mentioning 
that the conditions (dHJ) seem not to be comparable in any direct way with (fT5)) . and Nevertheless, our 
tools are totally different from those in [5] , [T^ and [15] . We prove our result by using the methods we established 
in [9] in which we proved the global well-posedness of the 31? viscous primitive equations. Finally, notice that 
the conditions ([1]) and (|14p are scaling invariant when the equation is considered on whole M^. Namely, 

2 3 
when — 'r — ~ 1, 




and 



a 13 
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hxWp dtj ^\^J^ \\p\r^ dtj , when - + ^ - 2, 

where ux = Xu{Xx, Xy, Xz, X^t) and px = X^u{Xx, Xy, Xz, X^t) which is also a solution to the Navier-Stokes 
equation if {u,p) is a solution to the Navier-Stokes equation in whole M.^. However, our condition is not scaling 
invariant. The reason may be because either our results are not optimal, or because our condition involves only 
one partial derivative of the pressure term and our result is limited to channel flows and not in the whole space 
M^, where the scaling argument is applied. 

The plan of this paper is as follows. In section[2l we reformulate the system (|2])- (|T0| . introduce our notations, 
and recall some well-known useful inequalities. Section [3] is the main section, and it is devoted for the regularity 
of solutions. 
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2. Preliminaries and Functional Setting 

In this section we introduce a new equivalent formulation of (|2|l- p0|) . Following the ideas introduced in [9] we 
integrate the equation (jj]) in the z direction and by ([7]) we obtain 



w{x,y,z,t) = - I Vh-v{x,y,£„t)d^, (17) 

and 







1 .1 

■ v{x,y,z,t)dz ■ v{x,y, z,t)dz = 0. (18) 

'0 Jo 

For every function 9{x, y, z), we denote by 

e{x,y)= f 0{x,y,z)dz. (19) 







and 



Following the geophysical terminology we denote by 



1 '■° 



vix,y)^- J ^v{x,y,OdL in M, (20) 
the barotropic mode. We will also denote by 

V = V — V, (21) 
the baroclinic mode, that is the fluctuation about the barotropic mode. Notice that 

9 = 0. (22) 

By substituting ([TT]) into we reach 

-^-uAhV~vv,, + {vVh)v~yj Vh-v{x,y,^,t)dn -£ + VhP = f, (23) 

and 

Vh-v = 0. (24) 

Furthermore, by taking the average of equations (|23p in the z direction, over the interval (0, 1) and using the 
boundary conditions ([7|), we reach 



-^~i^AhV + {vS7h)v~ U Vh-v{x,y,C,t)dn -^ + VhP^7- (25) 

As a result of (jlSp . (I22p . and integration by parts and using the boundary conditions (O, the nonlinear term in 
([25]) gives 



{v-Vh)v~{ I Vh-vix,y,^,t)d^]^ = iv-Vh)v + [iv-Vh)v + iVh-v)^. (26) 







By subtracting p5|) from ([23]) and using ([26|) we obtain 



— - vl^hV - vv^^ + {v ■ Vh)v - i Vh- v{x, y, ^, i)dM — 



■ V^)U + {v ■ \'h)v ~ [{v ■ Vh)v + {Vh-v)lI\ + VhP = /. (27) 
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In addition, v satisfies the boundary conditions: 

v{x + 1, y, z) = v{x, y + l,z)^ v{x, y, z), 

= 0. 



dv 

dz 



2 = 



dv 
oz 



(28) 
(29) 



z = l 



For convenience, we recall the following Gagiliardo-Nirenberg, Sobolev and Ladyzhenskaya inequalities in 
(see, e.g., [J, [H], [S], 124) 



"(M) < CaU\\%'lM)U\\Hi(^M) 



(30) 



for every (j> e H^{M), 2 < a < oo, and the following Gagiliardo-Nirenberg, Sobolev and Ladyzhenskaya inequali- 
ties in R3 (see, e.g., [J, [H], [H], [25) 



|L = (f2) < Ca||V'llL%2)IIV'll^/f(n) ' (31) 

for every ip e iJ^(ri), 2 < a < 6. Here Cq, denote constants which are scale invariant. Also, by ([50]) we get 



6-c. 3(q-2) 



<C|| 



Il°(a/) 



a/2||2//3 



,l"/2| 



l'/>l 



Q-2 



+ ll'/'l|L"(Af), 



(32) 



for every cj) € H^(M), and (3 > a. Also, we recall the integral version of Minkowsky inequality for the spaces, 
/? > 1. Let fii C M™^ and C M™^ be two measurable sets, where mi and m2 are two positive integers. Suppose 
that f{£,,ri) is measurable over 51i x Then, 

1 1//3 



< 



a2 \JOi 



1//3 



(33) 



3. Existence of the Strong Solution 

In this section we will prove the global existence of the strong solution to the system ([2|)- ([T0| under the 
assumption (|16p on . 

Theorem 1. Suppose that f,g £ H^{Q). For every {vo,wo) S V, and if p^ G L"([0, T], L^'?) with T > 0,a > 3 
and q> 1, then there is a unique strong solution {(v,w),p) of the system ©-(jl]) on [0,T]. 



Proof. In Remark [T] we described an argument that guarantees the existence of a Leray-Hopf weak solution and 
short time existence of the strong solutions. Suppose that {{v,w),p) is the strong solution with initial value 
ivo,wo) e V such that {v,w) e C{[0,T*),V) D L^{[0,T*), H^{n)), where [0,T*), forT* < T, is the maximal 
interval of existence. If T* = T, then there is nothing to prove. Next, we would like to show that certain norms 
of this strong solution remain finite for all the time, up to T*, provided the condition (|16p is valid. In this way 
we show that T* is, at least, equal to T. Namely, the strong solution {{v,w),p) exists on [0, T]. By the energy 
inequality ([72]) we have (see, for example, [12], [37], [40] for details) 

\\vmi + \\wmi<K,^, (34) 

and 
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where 



11/111 + ll.9lli 



1^2 A? 



+ \\vo\\l + \\wo\\ 



Ki2{t) = 



(ll/lii + llg|li)^ 



+ \\M\l + \\wo\\l 



(36) 
(37) 



3.1. ||v||r,3 < r < 4 estimates. Taking the inner product of the equation ([27|) with "^v in ^^(Sl), we get 
1 dlH; j ^|^^-,|2|^|,_2 + _ 2) I V;,|F| I' + \v,\^\vr^ + (r - 2) I 8^ l' l^^r~') ^^^ydz 



(u- V/i)?J 

n \ 



Qy 



-(v • Vh)v + (V;i • w) V + Vhp - /j • w dxdydz. 

By integration by parts and using the boundary conditions ([28]) and ([29|) we get 

i9u 



(w • V;i)t) 



• v{x, y, ^, j ^ ) • l-yT dxdydz = 0. 



(38) 



Moreover, ([24|) and the boundary condition ([29]) give us 



(u • Vft)w • |f I'' dxdydz = 0. 



(39) 



Thus, from dM])-® we have 
1 dll^JII 



+ / (iV/.iJI^I^fl'-^ + (r - 2) I Vh\v\ I' |?fr-2 + + (r - 2) I |" j^ir'^) dxdydz 



{v ■ V/Ju - (i; • \7h)v + (V;i • w) ^; + VhP - fj ■ I'l'T ^ dxdydz. 
Notice that by integration by parts and using boundary conditions (|^ and ([^ we have 
- y {v-Vh)v~ {v-\7h)v + {Vh-'^ V + V,,p] ■ l^r-^u dxdydz 

^ ~ ' -p(v,.(i^yr2iy)) 



k.j=l 



Observe that since p = 0, we have the Poincare inequahty 



dxdydz. 



\p\< \pz\ dz. 

J-h 



(40) 
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Therefore, from all the above and by Cauchy-Schwarz and Holder inequalities we obtain 

ld\\vrr 



r dt 

< C 

+C 

+c 

< c 



M 



+ y 1^ (|V^^jn^yr-2 + (r - 2) I Vu\v\ I' \vY-^ + \vA''\vY-^ + (r - 2) | d,\v\ |' l^l''-^) dxdydz 
dxdy 

dxdy 



1^1 / iv,.«i i«r-i dz 

Jq 

1 /.I 

2 J. / IV7 ~l l~|r-2 



M 



M 



\v\^ dz / |V,,w| \v\'-^ dz 
Jo 

Ip^I / |V,,t;| |i;r-2 ^^.^y _^ iijii^ ||„ 
Jo 

1/2 / ^1 X 1/2- 

2 l~|r-2 



ir-l 



N / \^hv\' \vr' dz 











i^r dz 



C [ f |wp dz ( f |V/,v|2 d^ ( f \vY'^ dz 

JM Jq \Jq J \Jq 



dxdy 



-C 



M 



\p.\dz ( C\Vhv\^ \vr^dz\ ' ( C W 















1/2- 



1/2' 



1/2 



< C||tJ|U4(M) ( j^^ \^hW \W~ dxdydz 



+C \ I [I dz ] dxdy 



+C\\\pMl^HM) [ 1^ I'^hvl' \vr^ dxdydz 



M \J0 



v\^ dz j dxdy 



dxdy 

dxdy + 11/11, 

\ 1/4 



V,iu|^ IwT"^ dxdydz 



v\^ ^ dz ) dxdy 



2(r+2) 



1/2 



M \J0 



1 \ 1 

r-2 



1/2?' 



dz da;dy 



where + l/g' = 1. By using Minkowsky inequality ((33)) . we get 

(iv (/ 1''''^'^^) '^^'^^j -"^Z (/ l^f'' dxdy^ ' dz. 
By virtue of (f32|) . we have 

/ dxdy < Co / |?^r dxdy [ \W^\Vhv\'^ dxdy + ( [ \v\'' dxdy 

JM JM JM \JM 

Thus, by Cauchy-Schwarz inequality we obtain 

/ 1 \ ^^'^ 1/2 
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Similarly, by ([33]) and ((32|) and (|30l) . we also obtain 

N 2/(r+2) 



M \J0 
1 



1 N ('■+2)/2 \ ^1 

2 > dxdy I < C 

"o V"'a/ 



\v\^ dz 



2/(r+2) 

dxdy ] dz 



< C 



< c 







L'-(M) 



\v\^-'\VhW dxdy] +\\v\\i'' 



'■(M) 



dz 



dz 



<C\\v\\l\\Vhv\[. 



r — 2 
7+2 



W~ l^^hvl dxdydz 



and 





r-l 


(HI 













(r+2)/(r-2) 



\ --+2 ^1 



\Ja/ 



wr-1V,,w|^ dxdy 



M 



2(r-2) 
'-('■ + 2) 



(r-2)/(r+2) 

\v\'^+'' dxdy ] dz 



|r-2 
L'-(Af) 



dz 



<C||??I1 



»■(■— 2) 



br"^|V/if;|^ dxdydz 



2(-— 2) 
r(r + 2) 



Thanks to (|5D|) and we conclude 



T--2 

1 \ 7^2 \ /.I 

r-2 



?rr-^ dz dxdy < C / / '^a;d2/ ) dz 

M \JQ J I Jq \Jm 



< C 



2(r-2) 



l|V/«t^llL2(M) + Ik 



|r-2 
Il2(m) 



dz < iiv/.^j|ir'" 



Il2(m)- 



Moreover, 



< C 



nr/q' 
\\L'-(M) 



\vr^\Vhv\'' dxdy 



M 



(r-2),'-r 



\l^{M) 



dz 



<C\\v\\;/'^' ( / \vr^\VHv\^ dxdydz 
In 



(r-2)g'-r 



(r-2) 
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Therefore, gT])-(|45l) and §0^ give 

d\\vrr 



dt 



+ v (|v„«n^;r-2^|V/,|iy||'H'-2 + |iy,|2|^;|'-2^|a,|«||'|^jr~2) dxdydz 



<C\\v\\l" \\Vuv\\l"\\v\\^'^ 



3/4 



\Vnv\' l^^r dxdydz ) + C\\v\\\'^ II V,TJ||^/^P1|; + ||/||. \\v\\l-' + C\\v\\ 



+C||^^|||||V,,?||2'- \\v\\r' ( / \v\'^-''\Vhv\'' dxdydz 

In 

r+2q' 



-C\\p.h, \\v\ 



V/ii;| dxdydz 



Ivq' 



\^hW \W~ dxdydz 



1/2 



By Young's and Cauchy-Schwarz inequalities we have 



dt 



(\Vhv\M'-^ + \^h\vfH 



-2 ,.|iJ j2|^|.-2^ 



\d,\vf \vY'-^^ dxdydz 



< C 



i + INl2llv.^||^ + ||^J|iniv,?||: 



l«li; + c|b.||^, + c||/||;. 



Thanks to Gronwall inequahty, we get 



II^^WIi; + ^/ / \\^y^?\vr^ + \'^AM^~^\ dxdydz <Kr, 



Jn 



where 



J^^ ^ ^CT+Ki^KMT)+K^ Ki2(T) 



1 + lko 115^1 (o)+ / Ib.(s)ll2,rfs + ||/||;r 



(46) 



(47) 



It is worth mentioning that the above estimate is also valid for 2 < r < 4. However, one need r > 3 in order to 
get the following estimate. 

3.2. estimates. Before we show the global bound, let us prove the following Lemma. 
Lemma 2. Let (p G H^{n),tp 6 L'^{n), and v,r be as in Theorem\^ Then, 

\v\ 101 IVI dxdydz < e (||V0||2 + + ||^||2) 

\^\v\\f-^ + Fll^ + (1 + \\v\\l) {\\v\\l + \\V^v\\l) 

for every e > 0. 
Proof. Notice that 

/ |d| \(t)\ 1-01 dxdydz < / (|w| + |w|) |0| |-0| dxdydz. 
Jn Jn 

By Holder inequality and ([50]) . we obtain 

l^l 101 1-01 dxdydz < \\v\\r ||0|| 2r 



(48) 
(49) 



<C||iy||. 110112- (||V0||2 + ||0.||2)N|0||2+N|r||0||2 

By Young's inequality, we reach 



\v\ 101 101 dxdydz < - (||V0||2 + 110 



z\\2 ■ 



+ cj\\v\\r' + \\v\\l] II0II 
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On the other hand, by applying the same method for provmg Proposition 2.2 in [8], we get 

\v\ \cj^\ 1^1 dxdydz < C (||iJ||^/'||V,tJ||^/' + \\vh) {mI^' ||V,<^||2 + Uh) Uh- 
Again, by Young's inequahty, we obtain 

\v\ 101 IVI dxdydz < I (||V,0||2 + ll^li^) + Q (1 + \\v\\l) {\\v\\l + \\Vhv\\1) UWl 
Therefore, El) holds. 



□ 

Next, let show the norm of the strong solution {v, w) is bounded. Taking the inner product of the equation 
([2]) with —AhV — Vzz and the equation ([3|) with — A^w — Wzz in , and using the fact that the Stokes operator 
is same as the Laplacian operator under periodic boundary conditions, we obtain 

l d{\\Vhv\\l + \\Vz\\l + liVfeWlli + 11^7,111) ||2^o||v7 ||2^|| ||2^||A ||2^o||v7 ||2^|| ||2N 

2^ + ^ y\\^hv\\2 + 2||V,,U^||2 + WVzzh + l|A/iW||2 + 2\\VhWz\\2 + W^zzh) 

= I [[v ■ Vh)v + wvz — /] • (A/iW + Vzz) dxdydz + [v ■ VhW + wwz — g] (Aft,w + Wzz) dxdydz. 



By integration by parts and using the boundary conditions ([7]) and (l8|), we obtain 

wVz ■ (AhV + Vzz) dxdydz 

= {{{VhW Vhv) ■ Vz + Wz\vz\'^) dxdydz 
In 

(VhWz ■ VhV + VhW ■ VhVz + Wzz Vz + Wz Vzz) ' V dxdydz, 



and also 



wWz {AhW + Wzz) dxdydz 



Vhw\'^Wz + wVhWz ■ VhW + -wl ) dxdydz 



2 

1 
2 



^ / (iVftWp + wl) Wz dxdydz 



(|V/iWp + wl) (V/i • v) dxdydz 
= 1 Vh {\Vhw\'^ + wl) ■ v dxdydz. 



Then, from the above, we get 
ld{\\Vhv\\l + \\vz\\l + \\Vhw\\l + \\wz\\l 



+ V (||A^w||2 + 2||V,,t;,|i2 + Wvzzf, + \\^hw\\ + 2\\VhWz\\ + \\wzz\\ 



2 dt 

= / {{'^ ■ V/i)w - /] • (A,iU + Vzz) dxdydz + / i^hWz ■ V hV + V hW ■ VhVz + Wzz Vz + Wz Vzz) ■ v dxdydz 
Jn 

+ [v ■ VhW 



j [v ■ VhW - g] {AhW + Wzz) dxdydz -\ ( "^h {\Vhw\^ + wl) ■ v dxdydz. 
Jn ^ Jn 
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By applying Lemma [2] with some small enough e and the Cauchy-Schwarz inequality, we obtain 

d{\\yhvM+\\Vz\\l + \\yhw\\l + \\w,\\l) ||2^ollT7 l|2_Lll l|2^||A l|2^0||V7 ||2^|| ||2N 

<c{\^ Uf"' + 11^^112 + l|V,.tJ||^) (||V,^;||2 + \\v,l\ + ||V,y;||2 + \\w^f^ + ||/||2 + ||5||2. 
Thanks to Gronwall inequality, we obtain 

\\'^Hvm\^hm?2 + \\^Hw{m^\\wMf2 

+v [ {\\Ahvis)\\l + 2\\yhv,{s)\\l + \\vAs)\\l + \\^hw{s)\\l + 2\\VhW,is)\\l + \\w,,{s)\\l) ds < K2, 



where 

K, = e^™^ (T+A-..(T)) + <<'-' T + + ||/||2 + 11^1,2] (50) 

Therefore, the norm of the solution remains bounded on the maximal interval of existence [0,T*). This 

completes the proof of theorem. □ 
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